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The present paper studies the influence of the structural modulation on the low energy physics 
of the Sr\4-xCaxCu24,OAi oxides, using ab-initio determination of the on-site and nearest neighbor 
effective parameters. The structural modulations appears to be the key degree of freedom, respon- 
sible for the low energy properties, such as the electron localization, the formation of dimers in the 
1 = compound or the anti-ferromagnetic order in the x = 13.6 compound. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.23.Ft 



The family of Srii-xCaxCu2iOii transition-metal ox- 
ides has attracted a lot of attention in the last decade. In- 
deed, the superconducting state observed in the x > 11.5 
compounds under high pressure , is supposed to be the 
realization of the remarkable theoretical prediction of su- 
perconductivity in two- legs doped spin ladders Q . More- 
over, this family of compounds exhibits a large diversity 
in electric and/or magnetic properties when chemical — 
isovalent substitution of Sr by Cb — and physical pres- 
sure are applied. For instance, under applied pressure, 
the compounds change from semi-conductor to conductor 
and finally to superconductor. 

This family of compounds possesses a complex layered 
structure of two alternating subsystems . The layers 
of the first subsystem are composed of weakly-coupled 
Ou02, spin-1/2 chains along the c direction. The spins, 
supported by 3c? orbitals of the Qi^+ ions, are coupled 
via two 90° Cu-O-Cu bonds. The layers of the second 
subsystem are composed of weakly-coupled, two-leg spin- 
1/2 ladders also along the c direction. The spins are 
strongly antiferromagnetically coupled on both legs and 
rungs, due to the 180° Cu-O-Cu bonds. The cell parame- 
ters of the two subsystems, in the direction of both chains 
and ladders, are incommensurate. The compounds have 
a pseudo-periodicity of 10 chain units for 7 ladder units. 

Electro-neutrality analysis shows that these systems 
are intrinsically doped with six holes by formula unit 
(f.u.). Similar to high-Tc superconductors, the holes 
are expected to be mainly supported by the oxygen 2p 
orbitals and to form Zhang-Rice singlets with the 
associated-copper hole. NEXAFS experiments 0| have 
later supported this assumption. A calculation of the 
Madelung potential on the concerned oxygen sites sug- 
gests that for the undoped compound, Sr 140^24041, the 
chains exhibit a larger electro-negativity than the lad- 
ders, resulting in a localization of all the holes on the 
former. Different experiments 01 , however suggest 
that about one hole per f.u. is located on the ladders. 
Under Ca substitution, the same experiments, as well as 
Cu NMR , show a transfer of part of the holes to the 
ladders. However the precise number of transferred holes 
is still under debate. X-ray data suggest a small hole 



transfer to the ladders (1.1 for x = 12) , while optical con- 
ductivity and ^^Cu NMR studies IS] show a larger hole 
transfer, of respectively 2.8 {x — 11) and 3.5 [x — 11.5). 

Let us first focus on the undoped system, Sri4Cu2404i. 
This is a semiconductor with a 0.18 eV gap. The spin 
ladders have a singlet ground state with a spin gap of 
about 35 — 47 meV 0, 0|- Surprisingly the spin 
chains also exhibit a si ngle t ground state with a spin 
gap of 11 - 12 meV [13 El [H E [11 . Since homoge- 
neous spin chains are known to be gap-less in the spin 
channel, the existence of a gap witnesses their strongly 
inhomogeneous character. In fact, the electronic struc- 
ture of the Sri4Cu2i04i chains is usually understood as 
formed by weakly interacting dimmers 
is now well established that these dimers are formed by 
second-neighbor spins separated by a Zhang- Rice singlet 
(ZRS), and order [ijj, ll2j llal according to a pseudo 
periodicity of 5 sites (one dimer followed by two ZRS). 

When Sr is substituted by Ca the system becomes more 
metallic. The doped compound with a; = 10 shows a gap 
of only 0.023ey 01 ■ This increased conductivity, as 
a function of the Ca doping, is usually understood as a 
consequence of the hole transfer from the chains to the 
ladders, in which the conduction is supposed to occur. 
However, a possible enhancement of the holes mobility 
within the chains is also evoked 0, 0|. Neutron scat- 
tering experiments jl8l | showed that the dimerization be- 
comes unstable with Ca substitution and disappears for 
X > 8, although the magnetic interactions within and 
between the dimers remain unchanged. In parallel, both 
ESR E3 ^'^d thermal expansion data 01 witness a pro- 
gressive disappearance of the charge order with increas- 
ing doping. Finally, at large doping {x > 11) and very 
low temperatures (< 2.5K) an anti- ferromagnetic phase 
is observed 0, 0| . 

An important aspect of these compounds, which is 
most of the time neglected, is the modulation of the 
two subsystems. Indeed, the mutual influence of the 
two subsystems results in a modulation of each of them 
with the periodicity of the other. These distortions are 
particularly large on the chains. Indeed, in the highly 
doped systems, the Cu-O-Cu angle varies between 89° 
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and 99°, while the Cu-0 distance varies with an am- 
plitude of 19%. The magnetic interactions being very 
sensitive to both bond angles and distances between the 
magnetic sites and the bridging ligands, one can expect 
that the modulations will be of importance for the low 
energy physics of the compounds. In the ladder subsys- 
tem, the structural distortions are of weaker amplitude 
since the alkaline-earth counter-ions are attached to it. 
In addition, the Cu-O-Cu angle varies around 9 = 180° 
and therefore super-exchange mechanism (that scales as 
cos^6' ~ 1 — 26^) should be dominant. Thus, while the 
effect of the structural modulations on the ladders may 
be of importance, it can be expected to be much weaker 
than on the chains sub-system. 

The aim of this paper is to study to which extent 
the modulations influence the chains electronic structure. 
For this purpose, we performed ab-initio calculations so 
that to accurately evaluate the influence of the modula- 
tions on the magnetic orbital energies (OE) and nearest 
neighbors (NN) interactions. We choose for these cal- 
culations the 5ri4Ctt2404i and 5rn.4Cai-i. 60^24041 com- 
pounds in their low temperature phase [20j. 

As the interactions between magnetic sites are essen- 
tially local, they can be accurately determined usin g an 
embedded fragment ab-initio spectroscopy method [21|. 
The fragment includes the magnetic centers, the bridging 
oxygens mediating the interactions, and their first coordi- 
nation shell. Short-range crystal effects are thus treated 
explicitly while the long-range crystal effects, such as the 
Madelung potential, are treated within an appropriate 
bath. The calculations have been performed using the 
DDCI method ^^l; ^ selected multi-reference single and 
double configuration interaction that properly treats (i) 
the strongly correlated character of the system, (ii) the 
mediation of the interaction via the bridging ligands and 
(iii) the screening effects on these processes 23j- The 
reference space has been chosen to be composed by the 
copper magnetic orbitals. The basis sets are of valence 
3C |2J| quality for Ou and 2C,+p for O. This method has 
proved to be very accurate in the determination of the 
the local effective interactions as well as local electronic 
structures. One can cite, for instance, the remarkable 
results obtained on the effective exchange and hopping 
determination in high- T,. p arent compounds [23, copper 
chains jl^ and ladder j^jj systems, where the computed 
values are within experimental accuracy, as well as on the 
charge ordering in the sodium vanadate low temperature 
phase 28]. 

Singlet-triplet excitation energies on embedded Cu20e 
fragments will thus yield the effective exchange integrals, 
while the first doublet-doublet excitation energies and 
associated wave functions yield both magnetic (hole) OE 
and hopping effective integrals. In order to study the 
influence of the incommensurability on these parameters, 
we performed calculations on 11 fragments associated to 
successive cells of the chain sub-system. 
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FIG. 1: Energy of the hole orbital as a function of its position 
along the chain, a) for a:: = 0, b) for x = 13.6. 



Let us first look at the hole OE (fig. 01. One sees im- 
mediately that their modulations are very large. Indeed, 
in the non-doped system, the OE vary within a range of 
1.2 eV. In the Ca doped system, the variation is even 
much larger and spans a 2.2 eV range. Let us notice, 
that in both systems, the OE variation is larger than 
the hopping and exchange energy scales (respectively of 
150 meV and 20 meV, see below). Thus, this is the OE 
that will dominate the low energy physics through the 
localization of the magnetic electrons (holes) on the low 
(high) energy sites. One should notice that the varia- 
tions of the OE are due to the crystal distortions. In- 
deed, when using the average crystal structure — where 
the subsystems distortions have been omitted — the OE 
exhibit only very small variations ^2^1 (< 30 meV for the 
undoped compound). 

Before studying the electrons (holes) localization along 
the chain, let us take a look at the variations of the NN 
hopping and exchange interactions. 
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FIG. 2: Exchange J and hopping t integrals as a function of 
the position of the fragment along the chain, a) J for a; = 0, 
b) J for X — 13.6, c) t for x = 0, d) t for x = 13.6. 

The effective exchange integrals, obtained for the 11 
fragments, are reported in fig.|21 As expected the NN ef- 
fective exchange for the x = compound is ferromagnetic 
and exhibits small modulations around an average value 
of 21.3 meV, the standard deviation being 2.5 meV. On 
the contrary, the x = 13.6 compound does not follow 
the expectations. Indeed, the effective exchange varies 
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greatly, going from ferromagnetic values (20 meV) up 
to anti- ferromagnetic interactions as large as —10 meV. 
These anti-ferromagnetic interactions are observed either 
for large Cu-O-Cu angles (~ 98°) or for large angles be- 
tween the two magnetic orbitals (> 10°) associated with 
short Cu-Cu distances. In both cases, these strong dis- 
tortions allow super-exchange mechanism to take place. 

The effective hopping integrals between NN copper 
atoms are expected to be quite small, since the nearly 
90° Cu-O-Cu angles forbid the through-bridge contribu- 
tion via the oxygen orbitals. However, this is not what is 
observed in our calculations (see fig-EJ since the hopping 
integrals present very large modulations and can reach 
values as large as 208 meV for x = and 266 meV for 
X = 13.6. These amplitudes are as large as 1/3 of the 
hopping observed in systems with 180° Cu-O-Cu angles, 
such as high-rp superconductors or CaO^ corner- 
sharing chains |23| . 

At this point it is clear that, the systems modulations, 
and the variations of these modulations according to the 
Ca doping, are crucial for the low energy properties of this 
family of iso-electronic compounds. We will now study 
whether they can explain the existence of the dimeriza- 
tion that is observed in the weakly doped compounds 
but not in the highly doped ones. For this purpose we 
need to extrapolate the values of the OE over the whole 
chain. It is usual to do so using a Bond Valence Sum 
(BVS) analysis, that analyzes the distances between the 
metal atom and its first coordination shell. However, for 
the 5ro.4Ca13.6O/24O41 system, the BVS analysis yields 
quite different results from the OE. A further analysis 
discloses that, in these systems, the Madelung potential 
on the magnetic centers, and thus the OE, is sensitive 
to the atomic displacements not only of the first coordi- 
nation shell, but up to the 8'^ shell of neighbors. This 
is for this reason that the BVS analysis fails to correctly 
reproduce the copper valence for this type of compounds. 
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FIG. 3: Orbital-energy differences as a function of the fourth 
coordinate r, a) for a: = 0, b) for x = 13.6. Solid line represent 
the fit. 



In order to extrapolate the OE over the whole chain, 
we used the crystallographic description of the incom- 
mensurate structure in a four dimensional space. Each 
subsystem is thus described by the three a, h and c usual 
spatial dimensions and a fourth coordinate r that has 
the periodicity of the other subsystem and describes the 



modulations |30| . We fitted the computed values using a 
Fourier analysis as a function of r (see fig. 13). One can 
note that the doped and undoped compounds present 
very different OE curves as a function of r. Let us notice 
on fig. the half periodicity of the magnetic cell for the 
SriiCu2AOAi cornpound, observed in neutrons scattering 
experiments [lllll8| . 

Figure0]reports the electrons (holes) localization along 
the chain as derived from the extrapolated on-site orbital 
energies. In the x — 13.6 CcLSG, Oil Fourier analysis has also 
been done on the effective NN exchange in order to pre- 
dict the sign of the NN interactions for the different posi- 
tions along the chain. The second neighbor exchange has 
been considered to be anti-ferromagnetic in agreement 
with experimental results. 

For the undoped Sri4Qi2404i system we studied the 
two types of filling considered in the literature, namely 
with all the holes on the chains (fig. EKl) and with one 
hole transfer per f.u. (fig. 0^2). When all the holes are 
on the chains, we retrieve a spin arrangement consistent 
with the experimental observations U 0, ^| , that is 
dimeric units formed of two second neighbors spins, sep- 
arated by two ZR singlets. These dimers are clustered 
by three or four units separated by a free spin, that is 
a spin with neither first nor second neighbor spin. We 
found 5 free spins for 100 sites (0.5 per f.u.) to be com- 
pared with the magnetic susceptibility measurements of 
0.55 free spins per f.u. The number of dimers is 17 
(1.7 per f.u.) to be compared with the magnetic suscep- 
tibility finding of 1.47 per f.u.. When one hole per f.u. is 
transfered to the ladders, the picture is totally modified. 
Indeed, the free spins totally vanish and first neighbor 
ferro-magnetically coupled spins appear. At the same 
time the number of (isolated) dimers is strongly reduced, 
with only 5 for 100 sites, that is 0.5 per f.u. From this 
analysis it is clear that if there are holes transferred to the 
ladders (at low temperature) in the undoped compound, 
this number is much smaller than 1 per f.u.. 

For the x — 13.6 doped system we studied three types 
of filling where n = 1,2, or 3 holes have been transferred 
to the ladders. For all fillings, we retrieve arrangements 
with nearest neighbors spins and no second neighbors 
dimeric units. The chains are essentially composed of 
low-spin clusters (either in a singlet or doublet state), 
with spin arrangements that present weak exchange frus- 
tration. Such fillings are thus expected to be specially 
stable, and should be put into perspective with the anti- 
ferromagnetic ordering seen in magnetic susceptibility 
and ESR measurements. For n = 1,2 one still observes 
a large number of free spins, while for n = 3 they have 
essentially disappeared, in agreement with magnetic sus- 
ceptibility experiments [l^ . 

In summary, we have studied the importance of the 
structural modulations on the low energy physics of the 
(^i4_2;Coa;Qi2404i family. Surprisingly these distortions 
are not simply responsible for parameters modulations 
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FIG. 4: Localization of the magnetic electrons along the chain, as a function of the number n of holes transferee! to the ladders, a) 
for X — 0, b) for x = 13.6. The number of sites considered in both cases corresponds to a better approximation of the compounds 
unit cells than the 10/7 ratio, that is c^h/cia = 3.9235/2.7268 ~ 100/69.5 for x = 0, and Cch/cia = 3.90136/2.73607 ~ 77/54 
for X = 13.6. Dots stand for ZR singlets, and the ellipsoids delimit spin clusters, which are magnetically very stable (with at 
most one second neighbor frustrated interaction). 



around their average value (except for the NN effective 
exchange on the undoped compound), but induce very 
large variations of the orbital energies, NN effective hop- 
ping and exchange. This is the variation of the OE 
(which spans a range of 1.2 for x = and 2.2 eV 
for X ~ 13.6) that is responsible for the low energy prop- 
erties of the compounds, through the localization of the 
magnetic electrons. It is in particular responsible of the 
formation of dimcrs in the undoped compound. In view 
of these results one can also suppose that the stabiliza- 
tion of part of the chains sites by larger structural mod- 
ulations, is responsible for the hole transfer toward the 
ladders in the doped compounds. In the x = 13.6 com- 
pound, the structural modulation is so large that it even 
reverse the sign of the NN effective exchange on part of 
bonds, lifting the exchange frustration that would arise 
for the large chain filling. In conclusion, one can say that 
the structural modulation is the key parameter, respon- 
sible for the large variation of the low energy properties 
in this family of compounds. 
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